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The influence of non-dimensional rotational velocity, flow Reynolds number and Prandtl number of the fluid on 


laminar forced convection from a rotating horizontal cylinder subject to constant heat flux boundary condition is 


numerically investigated. The numerical simulations have been conducted using commercial Computational Fluid 


Dynamics package CFX available in ANSYS Workbench 14. Results are presented for the non-dimensional rota- 


tional velocity a ranging from 0 to 4, flow Reynolds number from 25 to 40 and Prandtl number of the fluid from 


0.7 to 5.4. The rotational effects results in reduction in heat transfer compared to heat transfer from stationary 


heated cylinder due to thickening of boundary layer as consequence of the rotation of the cylinder. Heat transfer 


rate increases with increase in Prandtl number of the fluid. 


Keywords: Heated rotating cylinder, Non-dimensional rotational velocity, Flow Reynolds number, Prandtl num- 


ber 


Introduction 


The flow of fluids past heated stationary horizontal 
circular cylinder(s) is a challenging classical problem in 
fluid mechanics and heat transfer [1, 2]. Heat transfer 
from a stationary horizontal cylinder has been a subject 
of interest for many experimental and theoretical re- 
searches because of its numerous engineering applications. 

Dennis et al. [3] numerically investigated steady la- 
minar forced convection from a circular cylinder at low 
Reynolds numbers. The flow and energy equation were 
solved for Prandtl number up to 3.3x10 and Reynolds 
number up to 40 for the cylinder subjected to constant 
temperature. The predicted mean and local Nusselt num- 
ber was compared with the available experimental data of 
McAdams [4] and Van Der Hegge Zijen [5] and found to 
agree well. A numerical model was developed by Chun 
and Boehm [6] to study the forced convection heat trans- 
fer over a circular cylinder in cross flow. Two different 
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cases were studied viz. cylinder surface subjected to con- 
stant heat flux and constant wall temperature. The de- 
veloped numerical model were able to successfully pre- 
dict the heat transfer rate from the cylinder for Reynolds 
number up to Re = 2000, which were in good agreement 
with the laboratory experimental values. They also re- 
ported that uniform heat flux condition showed a higher 
value of heat transfer coefficient at the constant wall 
temperature of the cylinder for the same Reynolds num- 
ber. Experiments were performed by Sanitjai and 
Goldstein [7] to investigate the heat transfer by forced 
convection from a circular cylinder in cross flow for 
Reynolds number 2x10° to 9x10* and Prandtl number 
from 7 to 176. The cylinder was subjected to uniform 
heat flux boundary condition. The study reported that on 
the front part (0°< 8 < 85°) of the cylinder, local heat 
transfer increases with increasing Prandtl number. How- 
ever, on the rear part (85°< 6 < 180°), local heat transfer 
depends on flow characteristics near the surface. 
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Nomenclature 

C, specific heat of fluid (J/kg K) 

D diameter of the circular cylinder (m). 

h local convective heat transfer coefficient. 
(Wm°K `’). 

k thermal conductivity of fluid (W/m K) 

Lı upstream length from the inlet to the centre 
of the cylinder (m) 

Lə downstream length from the centre of the cy- 
linder to the outlet (m) 

Ly height of the computational domain (m) 

Nu local Nusselt number. 

Nüavg mean Nusselt number. 

P pressure (Pa) 

Pr Prandtl number. 

q" heat flux (Wm”). 

Re Reynolds number (U..D/v). 


Results of numerical investigation of momentum and 
heat transfer from cylinders in laminar cross-flow at 
10*< Re <200 was presented by Lange et al [8]. The fi- 
nite element code developed by Lange et al. [8] employ- 
ing multi grid and local grid refinement gave highly ac- 
curate results for drag coefficient, Strouhal and Nusselt 
numbers. Also for the Reynolds number considered in the 
study they were able to identify different flow regimes 
especially the critical Reynolds number where vortex 
shedding started was determined. Bharti et al. [9] ana- 
lyzed numerically heat transfer from an unconfined cir- 
cular cylinder in steady cross-flow regime. Numerical 
investigations were performed using a finite volume me- 
thod implemented on a Cartesian grid system in the range 
as 10<Re<45 and 0.7<Pr<400. A correlation for Nusselt 
number as a function of dimensionless variables was 
presented. The study revealed that for the same Reynolds 
number and Prandtl number the heat transfer coefficient 
from the cylinder surface is higher when it is subjected to 
uniform heat flux condition in comparison with uniform 
wall temperature condition. 

Recently forced convection heat transfer from rotating 
circular cylinders has gathered great importance for many 
researches due to its various applications which includes 
cylindrical cooling devices in plastics and glass indus- 
tries, contact cylinder driers in paper making, textile, and 
food processing industries, etc. The flow past such cy- 
linders is characterized by a series of wake interactions 
with increase in Reynolds number. These wakes interact 
with the hydrodynamic and thermal boundary layer over 
the heated rotating cylinder(s) altering the flow and 
thermal field. Badr et al. [10] numerically investigated 
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T temperature (K) 

Ta temperature of fluid at inlet (K). 

U Velocity vector. 

U» uniform velocity of fluid at inlet (ms"'). 

u x-component of velocity(ms”). 

v y component of velocity(ms`™”’). 

Greek letters 

a non-dimensional rotational velocity (Do /2U,,). 
v Kinematic viscosity (m7/s). 

P Density (kgm”). 

w Angular velocity (rads"). 

8 Angular displacement from the front stagnation 


point in the anti-clockwise direction (degrees). 


the laminar forced convection flow from a heated rotat- 
ing isothermal cylinder for Reynolds number up to 100. 
It was reported that for range of parameters considered in 
the study, the flow and thermal field were strongly influ- 
enced by the speed of rotation of the cylinder. Peramane 
and Sharma [11] studied numerically the flow and ther- 
mal field over a rotating cylinder for 20 < Re < 160 and a 
< 6. The study revealed that the mean Nusselt number 
and coefficient of drag decreases with increasing rotation 
speed. Numerical investigations were performed by 
Sharma et al. [12] to study the influence of Prandtl num- 
ber on flow and heat transfer from a heated rotating cy- 
linder subjected to constant wall temperature. They re- 
ported that the average Nusselt number increases with 
Prandtl number for affixed value of Reynolds number 
and non-dimensional rotational velocity. The above lite- 
rature review shows that many of the previous investiga- 
tors have analyzed the flow and heat transfer characteris- 
tics of heated stationary/rotating cylinders. In the above 
literature review it is also pointed out that the fluids in 
the vicinity of the heated cylinder can have a major in- 
fluence in the thermal energy transport process. It’s worth 
to mention here that the literature pertaining to the role of 
different fluids on the heat transfer from heated rotating 
cylinder is scarce. 

In the present work role of rotational rate of the heated 
cylinder and Prandtl number on the rate of heat transfer 
from a circular cylinder in the low Reynolds number re- 
gime in steady cross-flow is investigated. Detailed nu- 
merical results are presented for Reynolds number range 
25 < Re < 40, Prandtl number range 0.7 < Pr < 5.4 and 
non-dimensional rotational velocity range 0 <a <4. 
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Problem Statement and Mathematical Formula- 
tion 


2D, steady and incompressible flow of fluid of uni- 
form velocity U» and temperature T,, over a heated ro- 
tating circular cylinder of diameter D = 3cm is consi- 
dered (Fig.1). The cylinder is placed in an unconfined 
channel. This unconfined flow is simulated by placing 
the cylinder symmetrically between two walls which are 
subjected to slip boundary condition. The cylinder is 
placed at a distance of L;= 91.5cm from the inlet and at a 
distance L= 91.5cm from the outlet. The surface of the 
cylinder is subjected to constant heat flux of 100W/m’. 
Furthermore the thermo-physical properties of the fluid 
are assumed to be independent of temperature. The ob- 
jective of this work is to elucidate the role of Prandtl 
number, Reynolds number and non-dimensional rotation- 
al velocity on the forced convection from the heated ro- 
tating horizontal cylinder. The governing differential eq- 
uations pertaining to this problem is given below. 
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Fig. 1 Schematic diagram of the physical configuration 


Continuity Equation: 
V-U=0 
Momentum Equation: 


(U-V)U =- aww 
p 


Energy Equation: 
k 
(U-V)T =—— VT 
PC, 
Numerical Analysis 


The numerical analysis of the problem has been per- 
formed using commercial computational fluid dynamics 
package CFX available of ANSYS Workbench 14 [13]. 
The cylinder is modeled as a wall subjected to rotation 
and constant heat flux. The region between the cylinder 
and parallel plates is modeled as fluid. Simulations have 
been performed for different fluids of Prandtl number 0.7, 
0.9, 1.4 and 5.4. For each fluid simulations have been 
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performed for Reynolds number 25, 30, 35 and 40, a = 0, 
1, 2, 3 and 4. The maximum Reynolds number consi- 
dered in the study is less than 45. It is worth to mention 
here that as reported by Bharti et al. [9] for Re< 49 the 
flow can be treated as steady, incompressible and lami- 
nar. 


Boundary Conditions 


Free slip velocity condition is imposed on the outer 
walls. At the inlet uniform velocity is specified, wheras at 
the outlet outflow boundary condition is employed. The 
above boundary conditions is explained in detail as 
below. 

Inlet: 

Velocity- 

Outlet: 

The default outlet boundary condition in CFX, which 
assumes a zero diffusion flux for the flow variables. This 
physically means that the conditions of the outflow plane 
are extrapolated within the domain and have negligible 
impact on the upstream flow conditions. This is similar to 
Neumann boundary condition [14] 

Cylinder surface: 

No-slip boundary condition. Velocity- j=constant. 
Heat flux = 100W/m”. 

Outer wall: 

a =0,v=0 and i =0 
oy oy 

The local and the average Nusselt number are eva- 
luated by the following expressions. 


u=U,,, v=0. Temperature- T= To, 


hD 
Nu; = pa 
1 2m 
Nuayg = sale Nutigcai4 O 


Computational Domain and Grid 


A two dimensional computational domain has been 
created for numerical computations. The computational 
domain is discretized using structured non-uniform hex- 
ahedral cells. At the solid liquid interface and also at the 
region downstream of the cylinder the gradient of tem- 
perature and velocity are high and therefore to resolve 
these gradients fine grids have been employed in these 
regions. The grid is generated using the commercial grid 
generator ICEM CFD available in ANSYS Workbench 
14. In order to obtain an optimum grid for numerical 
computations a grid independence study was conducted 
and the results are presented in Table 1. It can be ob- 
served from Table 1 that a structured, non-uniform grid 
structure with total 1,12,800 elementsis found to be op- 
timum for numerical simulations. The final computation- 
al domain along with the grid used for the analysis is 
shown in Fig. 2. 
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Table 1 Grid Sensitivity Study (q” = 100 W/m’, a = 3, Re = 40, and Pr = 0.7) 


Numerically predicted surface temperature at different angular positions(0) 


No of cells 
Tı T2 T; T4 Ts To T: Ts To Tio Tu Ti 

28,200 337.6 334.4 331.8 330.6 330.6 331.4 332.6 334.0 335.5 337.1 338.5 338.9 

56,400 336.8 334.2 331.6 330.5 330.4 331.0 332.1 333.4 334.7 336.2 337.4 337.8 

1,12,800 338.1 334.9 332.2 331.0 330.9 331.7 333 334.5 336.0 337.6 339.0 339.4 

2,25,600 338.4 334.8 332.4 331,1 331 332.1 333.4 334.6 336.1 338.5 339.2 339.6 
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Close up view in the vicinity of heated cylinder 


Fig. 2 Computational domain along with grid used for 2D 


numerical simulation 


Numerical Solution Procedure 


The governing equations of flow and heat transfer are 
solved by using a finite volume solver CFX. The semi 
implicit method for the pressure linked equations 
(SIMPLE) has been used as the pressure velocity coupl- 
ing scheme. The upwind scheme is employed to discret- 
ize the convective and diffusive terms. Convergences of 
the solution are checked by monitoring the residues of 
discretized governing equations. A convergence criterion 
of 10° is used for continuity, momentum and energy eq- 
uations. 


Results and Discussions 


The numerical methodology employed in the present 
study has been validated with the previous numerical 
work reported in literature (Bharati and Chhabra) [9], to 
investigate heat transfer from a stationary heated cylinder 
in cross-flow. The comparison of results of Local Nusselt 
number variations shown in Fig. 3 is of good agreement 
with each other. 

To get more insight into the quantitative measure of 
heat transfer, a representative comparison of the thermal 
field is made and the result is shown in Table 2. The re- 
sults presented in Table 2 indicate that the average sur- 
face temperature is found to agree well with the work of 
Bharati and Chhabra [9]. 


circumference of the cylinder for different angular po- 
sitions (0). (Re = 20, q” = 100 W/m”, Pr = 0.7, a = 0) 


Table 2 Comparison of average surface temperature with 
literature value (Re = 40, q” = 100 W/m”, Pr = 0.7, a = 0) 


Numerically predicted 


No of elements 
surface temperature 


Bharti and Chhabra [9] 


28,200 - 331.3 
56,400 329.9 330.9 
1,12,800 - 332.4 


The good agreement between the results of present 
work and that of [9] is indicative of the efficacy of the 
numerical methodology employed here for the numerical 
analysis of heat transfer problems of similar nature. The 
numerical methodology is further extended to study the 
influence of rotation of the heated cylinder. Numerical 
simulations were performed for different non-dimensional 
rotational velocity (a = 0, 1, 2,3,4) at Re = 40, 35, 30 and 25. 

The simulation exercises were also performed to study 
the effect of Prandtl number (Pr = 0.7, 0.9, 1.4 and 5.4) 
that corresponds to fluids air, ammonia, methane and 
water respectively, while keeping the heat flux as 100 
W/m’. It is important to mention here that the rotation 
direction is taken in the anticlockwise sense. 


Effect of rotation on heat transfer 


The variation of local Nusselt number along the cir- 
cumference of the heated cylinder computed at different 
angular positions for various non-dimensional rotational 
speed corresponding to Pr = 0.7 and Re = 40 is depicted 
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in Fig. 4. Figure illustrates that as rotational velocity 
changes the position of maximum heat transfer shifts in 
the direction of rotation. As rotational speed increases the 
flow becomes asymmetric under the influence of the ro- 
tation of the cylinder due to the interaction between the 
flow generated by the rotation of the cylinder and the free 
stream flow. The streamlines and isotherms shown in 
Figs. 5 (a) and (b), respectively, clearly indicate the 
asymmetric nature of flow created by the rotation. The 
impact of rotation of cylinder results in thickening of 
boundary layers, therefore the transport of momentum 
and energy by diffusion is expected to be reduced. Con- 
sequently the rotation effects causes reduction in heat 
transfer when compared to heat transfer from the cylinder 
when it is in stationary condition. 

Interestingly, as seen in velocity streamlines depicted 
in Fig 5(a) the stagnation point shifts from the front end 
(equator) of the cylinder to the pole (top) as rotation 
speed increases. This causes the isotherm to lose its 
symmetrical nature and shift in the direction of rotation 
of the cylinder that is towards the bottom side of the 
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Fig. 4 Local Nusselt number (Nur) variation along the cir- 


cumference of the cylinder for different angular posi- 
tions (0) for air. (Re = 40, q” = 100 W/m’, Pr = 0.7) 


(b) 


Fig. 5 Streamline contour and isotherms (Re = 40, q” = 100 
W/m’, Pr = 0.7) 


cylinder. Fig. 5(b) conveys that increasing rotational 
speed causes the location of denser regions of isotherms 
to shift in the anticlockwise direction. Therefore, maxi- 
mum value of Nusselt number also changes from front 
stagnation point towards the bottom of the cylinder in the 
anticlockwise direction. 


Effect of flow Reynolds number on heat Transfer 


The influence of Reynolds number on heat transfer is 
studied by varying the free-stream velocity, while keep- 
ing non-dimensional rotational velocity and Prandtl 
number of the fluid at fixed values. The result presented 
in Fig. 6 shows that for a given dimensionless rotational 
velocity and Prandtl number, the average Nusselt number 
increases with flow Reynolds number. For a given di- 
mensionless rotational speed, as the flow Reynolds 
number increases the inertia force increases. Obviously, 
higher inertial force of the incoming fluid counteracts the 
flow field generated by the rotation thereby thinning the 
boundary layer. As a result the heat transfer increases due 
to increase in flow Reynolds number at a given rotation 
speed. Also it is important to note that increase in inertia 
causes the streamlines to bent more towards the heated 
cylinder (Fig. 7) resulting in higher heat transfer from the 
cylinder. Revisiting Fig. 6, it can be inferred that for a 
given flow Reynolds number the average Nusselt number 
decreases with rotational speed, the reason for this was 
explained in section 3.1. 


ed 
0 


270° 
—e— Re=40 r 
eO Re=35 gol i \igoo 
----- Re=30 wy 
= a Re=25 Pri 


Average Nusselt number, Nua 


Non-dimensional rotational velocity, a 


Fig. 6 Average surface Nusselt number (Nuavg) at different 
rotational velocity and Reynolds number. (Pr = 0.7, q” 
= 100 W/m’) 


(b) Re =25 


(a) Re =40 


Fig. 7 Streamline contour (a = 4, q” = 100 W/m”, Pr = 0.7) 
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Effect of Prandtl number on heat transfer 


In order to study heat transfer from a rotating cylinder 
in various fluids, numerical simulations were carried out 
with fluids having Prandtl number of different values 
such as 0.7, 0.9, 1.4 and 5.4. The results of study have 
shown that as Prandtl number increases the average Nus- 
selt number increases as presented in Table. 3. As we know, 
as thermal boundary layer thickness decreases with in- 
crease in Prandtl number the heat transfer rate is expected 
to increase with increase in Prandtl number. Moreover, 
higher momentum diffusivity of the fluid enhances the 
momentum transport by diffusion in the boundary layer. 
These supplements to the higher heat transfer rates. 

A comparison of the variation in local heat transfer 
coefficient for different angular positions along the cir- 
cumference of the cylinder plotted in Fig. 4 and Fig. 8 
shows marked difference in behavior. When the heated 
cylinder is stationary and placed in air the heat transfer 
decreases progressively from the leading edge to angular 
position at 180°, thereafter heat transfer increases. The 
reason for this characterization is evident from the plot of 
isotherms shown in Fig. 9. However, for heated statio- 
nary cylinder placed in water the heat transfer decreases 
to the point at 150°, followed by a small increase at 180°, 
thereafter the pattern repeat till the back leading edge is 
reached. This behavior is due to the shift in the isotherm 
pattern depicted in Fig. 9. Similar results were also ob- 
tained for other fluids but the deviations tend to decrease 
for smaller Prandtl number. 


Table 3 Comparison of Average surface Nusselt number (Nuavg) 
for different Prandtl number (Re = 40, q” = 100 W/m”, a = 3) 


Prandtl number (Pr) Average surface Nusselt number(Nuavg) 
0.7 3.23 
0.9 3.46 
1.4 3.66 
5.4 6.25 
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Fig. 8 Local Nusselt number (Nur) variation along the cir- 


cumference of the cylinder for different angular posi- 
tions (0). (Re = 40, q” = 100 W/m”, Pr = 5.4) 
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Fig.9 Isotherms. (Re = 40, q” = 100 W/m’, a = 0) 
Conclusion 


Numerical simulations have been performed to study 
the effect of non-dimensional rotational velocity, flow 
Reynolds number and Prandtl number on the flow and 
heat transfer from a heated rotating cylinder subjected to 
constant heat flux. The salient conclusions from the study 
are 
e — As rotational speed increases the location of maxi- 

mum heat transfer shift towards the bottom of the 
cylinder along anticlockwise direction due to the 
redistribution of flow and thermal field. 

e For a given flow Reynolds number, as rotational 
speed increases the thickness of the boundary layer 
developed due to the rotation of heated cylinder in- 
creases, thereby, deteriorating heat transfer per- 
formance. 

e While keeping the rotational speed constant, in- 
crease in flow Reynolds number causes stream lines 
to shift more towards the heated cylinder resulting 
in higher heat transfer rate. 

e — For low and high Prandtl number fluids considered 
in the study, heat transfer characteristics of a sta- 
tionary heat cylinder shows distinct features. 

e For a given flow Reynolds number, the variation of 
local Nusselt number is found to be quashing pro- 
gressively as non-dimensional rotational velocity 
increases. 
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